Three-body forces and shell structure in calcium isotopes 
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Understanding and predicting the formation of shell structure from nuclear forces is a central 
challenge for nuclear physics. While the magic numbers N = 2, 8, 20 are generally well understood, 
N — 28 is the first standard magic number that is not reproduced in microscopic theories with two- 
nucleon forces. In this Letter, we show that three-nucleon forces give rise to repulsive interactions 
between two valence neutrons that are key to explain 48 Ca as a magic nucleus, with a high 2 + 
excitation energy and a concentrated magnetic dipole transition strength. The repulsive three- 
nucleon mechanism improves the agreement with experimental binding energies. 



PACS numbers: 2l.lQ.-k, 21.30.-x, 21.60.Cs, 27.40.+Z 

In nuclei certain configurations of protons and neu- 
trons (nucleons) are observed to be particularly well- 
bound. These closed-shell or "magic" nuclei form the 
basis of the nuclear shell model (exact diagonalizations 
in spaces based on the observed shell structure [l| ) , which 
provides a key computational method in nuclear physics. 
Exploring the formation of shell structure and how these 
magic configurations evolve with nucleon number to- 
wards the drip lines is a frontier in the physics of nu- 
clei, and a microscopic understanding from nuclear forces 
presents a major challenge for theory. 

The theoretical shortcomings in predicting shell struc- 
ture are particularly evident in the calcium isotopes. 
While microscopic calculations with well-established two- 
nucleon (NN) forces reproduce the standard magic num- 
bers N = 2, 8, 20, one of the most striking failures is that 
they do not predict Ca as a doubly- magic nucleus when 
neutrons are added to 40 Ca [J, making N = 28 the 
first standard magic number not reproduced in micro- 
scopic NN theories. As a result, phenomenological forces 
have been adjusted to yield a doubly-magic 48 Ca [1,0]' 
and it has been argued these phenomenological adjust- 
ments may be largely due to neglected three-nucleon (3N) 
forces [a]. Recently, we have shown that 3N forces play 
a decisive role for the oxygen anomaly and can explain 
why 24 is the heaviest oxygen isotope [fj|. In this Letter, 
we present the first study of the impact of 3N forces on 
medium-mass nuclei. Our results demonstrate that one- 
and two-body contributions from 3N forces to valence 
neutrons, as well as extended valence spaces, are essen- 
tial to understand shell structure in the calcium isotopes 
and N — 28 as a magic number based on nuclear forces. 

Three-nucleon forces were introduced in the pioneering 
work of Fujita and Miyazawa (FM) Q and arise because 
nucleons are finite-mass composite particles that can also 



be excited by interacting with other particles. The long- 
range part of 3N forces is dominated by the FM 3N mech- 
anism, where one nucleon virtually excites a second nu- 
cleon to the A(1232MeV) resonance, which is de-excited 
by interacting with a third nucleon. Additional long- 
range and shorter-range 3N interactions are included nat- 
urally in chiral effective field theory (EFT) [8|, which 
provides a systematic expansion for nuclear forces. The 
importance of chiral 3N forces has been well established 
in light nuclei with A = N + Z < 12 

We derive the interactions among valence neutrons fol- 
lowing two approaches. First, we use low-momentum 
interactions V\ ow k with smooth cutoffs 11011 obtained by 
evolving a chiral N 3 LO NN potential [llj to lower res- 
olution with A = 2.0 fm -1 . The two-body interactions 
in the pf and pfgg/2 shell are calculated to third or- 
der in many-body perturbation theory (MBPT) follow- 
ing Refs. [H, HH in a space of 13 major shells. We use 
a harmonic oscillator basis with hui = 11.48 MeV, ap- 
propriate for the calcium isotopes. Our results are con- 
verged for intermediate-state excitations to 18hu. Sec- 
ond, we take a standard G-matrix that has been used 
as a starting point in many nuclear-structure calcula- 
tions [HI . This G matrix is based on the Bonn C NN po- 
tential and includes many-body contributions to third or- 
der, but with intermediate-state excitations to 2huj (and 
calculated for haj = 10 MeV). Although the G-matrix 
calculation should be improved by extending the inter- 
mediate states, this two-body interaction has been used 
as standard starting point for shell-model studies, and 
gives us a baseline to investigate changes due to three- 
body forces. A detailed study of the cutoff dependence, 
which provides a measure of the theoretical uncertainty, 
is left to future work. 

To study the validity of the MBPT approach, we 
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FIG. 1. Neutron SPEs relative to the 40 Ca energy as a function of mass number A. (a) SPEs obtained from phenomenological 
models KB3G || and GXPF1 [J. (b) NN-only theory: SPEs calculated from a G- matrix and low-momentum interactions 
Viowfc- (c) SPEs including contributions from 3N forces due to A excitations, 3N(A), and chiral EFT 3N interactions at N 2 LO, 
3N(N 2 LO) 0. The results in panels (b) and (c) start from the empirical GXPF1 SPEs in 41 Ca. (d) SPEs in the pfg 9/2 
shell calculated from Viowfc and 3N(N 2 LO) forces, where the starting SPEs in 41 Ca are calculated consistently in many-body 
perturbation theory [MBPT] and include one-body 3N contributions. 



have carried out coupled-cluster (CC) calculations for the 
ground-state energies of 4o,48,52,54 Ca at the A-CCSD(T) 
level [Hj], based on the same VWfc interaction and ba- 
sis space. We have verified that the CC energies are 
converged in these spaces. Using particle-attached (to 
40 Ca) CC energies as single-particle energies (SPEs), the 
MBPT results agree within a few percent with CC theory 
(relative to 40 Ca): -159.3 MeV (MBPT) vs. -155.0 MeV 
(CC) for 48 Ca; -230.7 MeV vs. -235.9 MeV for 52 Ca; and 
-259.1 MeV vs. -268.6 MeV for 54 Ca. This shows that 
MBPT can be comparable to CC theory for V\ ow k inter- 
actions, but also highlights the important role of SPEs. 

To understand shell structure in the calcium isotopes, 
we show in Fig. Q] the change of the SPEs of the neu- 
tron orbitals, with standard quantum numbers lj, in 
the pf and pfgg/2 shell as neutrons are added to 40 Ca. 
Fig. [2 ( a ) gives the evolution obtained from phenomeno- 
logical models: GXPF1 jij], a quasi-global fit of two-body 
matrix elements (starting from a G-matrix) and of initial 
SPEs in 41 Ca to experimental data; and KB3G Q, which 
modifies the monopole part of a standard G matrix. De- 
spite different initial SPEs, as neutrons fill the f-j/2 shell, 
the repulsive interaction between the f^/2 and p 3 / 2 neu- 
trons causes a significant gap to develop at V = 28, in- 
dicative of a shell closure. 

In contrast, both NN-only theories in Fig.[T](b) exhibit 
minimal repulsion between these orbitals and the gap re- 
mains largely unchanged from 40 Ca. This is even more 
evident when the 41 Ca SPEs are calculated microscop- 
ically at the same third-order level: the NN-only SPEs 
(based on Vi ow fc) are too bound and both / 7 / 2 and p 3 / 2 
orbitals are at -10.8 MeV. Starting from the GXPF1 
SPEs in 41 Ca, the NN-only results in Fig. [1] (b) depend 
only weakly on the approach or NN forces used, except for 
differences for N > 28 in the interactions among the p or- 
bitals, mostly £> 3 /2 — P3/2 (due to more attractive second- 
order core-polarization and third-order TDA/RPA con- 
tributions with Vlow/c, A = 2.0fm _1 , compared to the 
G-matrix) . This uncertainty remains when 3N forces are 



added [see Fig. [2(c)]. Beyond 48 Ca, both NN-only the- 
ories predict a shell closure at N — 34, which is a major 
disagreement between the phenomenological models. 

The dominant differences between phenomenological 
forces and NN-only theory can be traced to their two- 
body monopole components, which determine the aver- 
age interaction between orbitals @, [l5[ • In an operator 
expansion, the monopole interaction corresponds to the 
term involving number operators, so that differences are 
enhanced with TV, and the SPE of orbital lj is effectively 
shifted by the monopole interaction multiplied by the 
number of neutrons in orbital . The interplay of this 
two-body effect with the initial SPEs largely determines 
the formation of shell structure. 

Next, we include 3N forces among two valence neu- 
trons and one nucleon in the core. In Ref. 0], we have 
shown that these configurations give rise to repulsive 
monopole interactions among excess neutrons. This cor- 
responds to the normal-ordered two-body parts of 3N 
forces, which was found to dominate in coupled-cluster 
calculations [lij over residual 3N forces (the latter should 
be weaker because of phase space arguments [13] )■ For 
the G-matrix approach, we include FM 3N forces due to 
A excitations, 3N(A), where the parameters are fixed 



by standard pion-N-A couplings |18|. For chiral low- 
momentum interactions, we take into account chiral 3N 
forces at N 2 LO [l9(. These include long-range two-pion- 
cxchange parts Cj (due to A and other excitations), plus 
shorter-range one-pion exchange cd and 3N contact ce 
interactions that have been fit to the 3 H binding energy 
and the 4 He matter radius [licj . For Vi ow fe, we also con- 
sider the one-A excitation 3N force that corresponds to 
particular values for the two-pion-exchange part c; and 
c D = c E = [Vi owfc +3N(A)] [8]. In the Vi ow fe+3N calcu- 
lations, we also scale all matrix elements by hui ~ A^ 1 ^ 3 . 

For all results, full 3N multipolc contributions are in- 
cluded to first order 2l|, although only the monopole 
part is responsible for the SPE evolution in Fig. [1] (c). 
Here we see in both microscopic approaches, 3N forces 
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FIG. 2. First 2 + excitation energies in the even calcium 
isotopes compared with experiment, (a) Energies obtained 
from phenomenological models KB3G || and GXPF1 [J. 
(b) NN-only theory: energies based on a G-matrix and low- 
momentum interactions Vj ow fc with empirical GXPF1 SPEs in 
41 Ca, as well as with KB3G values [SPEJCB3G]. (c) Including 
contributions from 3N forces due to A excitations, 3N(A), and 
chiral EFT 3N interactions at N 2 LO, 3N(N 2 LO) (d) En- 
ergies from Viowfe and 3N(N 2 LO) forces in the pfgg/2 shell 
with empirical GXPF1 SPEs and g g/2 at -lMeV in 41 Ca, as 
well as with SPEs in 41 Ca calculated consistently in MBPT. 

provide repulsive shifts of all single-particle levels, chang- 
ing the binding energies as shown in Fig. |4j In addition, 
the repulsion between the fa/2 and P3/2 orbitals leads 
to an increased separation at N = 28, similarly to the 
phenomenological forces. Moreover, the gap at N = 32 
is increased due to the repulsive p 3 / 2 — P1/2 interaction, 
while N = 34 remains approximately the same. 

We take into account many-body correlations by di- 
agonalization in the valence space and plot the first 2+ 
energy of the even calcium isotopes in Fig. [5] The excita- 
tion energies of the phenomenological models in Fig. [2] (a) 
show the fit to the high 2 + energy in 48 Ca and hence the 
doubly-magic nature, and highlight the difference in the 
prediction of N = 34 as a shell closure. In contrast, 48 Ca 
is not reproduced in any calculation based on NN forces 
in Fig. [5](b), regardless of starting SPEs, or whether we 
include the g g / 2 orbit in Fig. [5] (d) . 

With 3N forces in Fig. H (c), the 2+ energy in 48 Ca 
is uniformly improved. The pf shell predictions with 
NN+3N forces are similar with initial GXPF1 SPEs, but 
still below the experimental value. With KB3G SPEs 
in 41 Ca, the 2 + energy is significantly lower due to the 
smaller initial fa/2 — P3/2 spacing. When the g g / 2 orbit is 
included in Fig.[2](d), the 2 + energy is obtained very close 
to experiment. In addition, we find that all microscopic 
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FIG. 3. Magnetic dipole transition rates from the ground 
state to 1 + excited states in 48 Ca compared with experi- 
ment [2J. The B(M1) values are calculated in the pf and 
pfgg/2 shell in panels (a) and (b), respectively, based on NN- 
only interactions and including 3N forces (spin g factors are 
quenched by 0.75). The results are labeled as in Fig. [2] 

NN-only and NN+3N results at this level yield a high 2+ 
energy in 54 Ca, and hence a shell closure at N=34 (as 
suggested in Ref. [22|). The similarities of V\ ow fc+3N(A) 
and +3N(N 2 LO) in Fig. [2] (c) demonstrate that the con- 
figurations composed of valence neutrons probe mainly 
the long-range parts of 3N forces. 

To remove the uncertainty in the initial SPEs, we cal- 
culate the SPEs in 41 Ca by solving the Dyson equation, 
consistently including one-body contributions to third or- 
der in MBPT in the same space as the two-body inter- 
actions, and chiral 3N forces between one valence neu- 
tron and two core nucleons to first order. In contrast 
to the failure with NN-only forces, we find in Fig. [1] (d) 
the pf shell SPEs are generally similar to the empiri- 
cal ones, and we find the gg/2 to initially lie between 
the P1/2 and f$/ 2 orbitals. Our results based on MBPT 
SPEs and consistent two- valence-neutron interactions are 
shown in Fig. [2] (d) . The agreement with experiment is 
very promising for a parameter-free calculation based on 
NN and 3N forces. Furthermore, the high 2 + in 48 Ca, 
despite a relatively small fa/ 2 — P3/2 g a Pj reflects the 
possible importance of correlations beyond the p/-shcll 
(in the context of SPEs, see also Ref. [23j). Another chal- 
lenge for microscopic theories is the prediction of the first 
excited (l/2~) state in 49 Ca, which indicates the size of 
the p 3/2 -P1/2 gap at N = 28. For VW fe +3N(N 2 LO) in 
both the pf and p/39/2 shell, this energy is « 1.0 MeV 
compared to the experimental value 2.02 MeV, while the 
MBPT results yield 1.8 MeV. 

We further examine the closed-shell nature of 48 Ca 
in Fig. [21 which shows the magnetic dipole transition 
rates B(M\) from the + ground state to 1 + excited 
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FIG. 4. Ground-state energies of calcium isotopes relative 
to 40 Ca compared with experiment and extrapolated energies 
from the AME2003 atomic mass evaluation [2q] . The panels 
and results are labeled as in Fig. [2] 



forces for binding energies and evolution of shell struc- 
ture in medium-mass nuclei, thus linking the 3N forces 
frontier to the experimental frontier for neutron-rich nu- 
clei. Our results show that 3N forces and an extended va- 
lence space are key to explain the N = 28 magic number, 
leading to a high 2 + excitation energy and a concentrated 
magnetic dipole transition strength in 48 Ca. It is intrigu- 
ing and promising that the parameter-free MBPT results 
in the extended valence space reproduce experiment best. 
Future work will include a detailed comparison to empiri- 
cally adjusted interactions, where the pfgg/2 interactions 
can also be transformed into p/-shell-only interactions by 
an Okubo transformation. 

This work was supported by the US DOE Grant DE- 
FC02-07ER41457 (UNEDF SciDAC Collaboration) and 
DE-FG02-06ER41407 (JUSTIPEN), by grants-in-aid for 
Scientific Research (A) 20244022 and (C) 22540290, the 
JSPS Core-to-Core program EFES, and the Alliance 
Program of the Helmholtz Association (HA216/EMMI). 
Part of the numerical calculations have been performed 
on Kraken at NICS, UT/ORNL, and at the JSC, Jiilich. 



states, where the experimental concentration of strength 
indicates a single-particle transition 24|. With NN-only 
forces, there is a significant fragmentation of strength, 
and the energy of the dominant transition is below the 
observed value. When 3N forces are included, the peak 
energies are pushed up for all 3N(N 2 LO) cases. More- 
over, the MBPT results predict a clear concentration, in 
very good agreement with experiment. 

Finally, we turn to the ground-state energies in Fig. [H 
which have been measured to 52 Ca [25| and are known 
to exist to 58 Ca [26| , With phenomenological models, 
the ground-state energies decrease to N=34, then the be- 
havior flattens to N=40 due to the weakly-bound / 5 / 2 
orbital. With NN-only forces in Fig. 2] (b) [as expected 
from Fig. [1] (b)], all neutron-rich calcium isotopes are 
overbound. In Fig.|4](c) and (d) the repulsion due to 3N 
forces leads to less bound ground-state energies, and all 
calculations with 3N forces exhibit good agreement with 
experiment [in Fig. H (c) the Mow fc+3N(N 2 LO) (KB3G 
SPE) results would lie on those of Mowfc+3N(A) (GXPF1 
SPE)]. The repulsive 3N mechanism, discovered for the 
oxygen anomaly [6j, is therefore robust and general for 
neutron-rich nuclei. In our best calculation with MBPT 
SPEs in the pfgg/2 shell, the ground-state energies are 
modestly more bound. Our results with 3N(N 2 LO) sug- 
gest a drip line around 60 Ca, which is close to the ex- 
perimental frontier [26j |. As the predicted energies can 
significantly flatten from N — 34 — 40, as is the case in 
our best MBPT calculation, the inclusion of continuum 
effects will be very important. 

We have presented the first study of the role of 3N 
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